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Geomagnetic Probing of Core Geodynamics 
and Earth’s Magnetic Shield 

As Earth’s main magnetic field weakens, our magnetic shield against the onslaught of the 
solar wind thins. And the field strength needed to fend off battering by solar coronal 
mass ejections is decreasing, just when the delicate complexity of modem, vulnerable, 
electro-technological systems is increasing at an unprecedented rate. Recently, a working 
group of distinguished scientist from across the nation has asked NASA’s Solid Earth and 
Natural Hazards program a key question: What are the dynamics of Earth’s magnetic field 
and its interactions with the Earth system? 

Some answers to part of this question come from diverse studies of Earth’s core 
geodynamo, a dynamical system as complex as weather and climate, yet buried deep 
under 1800 miles of rocky mantle. Electric current flowing in Earth’s metallic core is the 
accepted main source of the main field, but it is weakening faster than expected from the 
small electrical resistivity of the metal. According to one theory, as Earth cools over 
geologic time, the solid iron inner core freezes from the molten outer core. This liberates 
latent heat and buoyant slag which rise, stirring the outer core to motion at speeds of 
several miles per year. But whatever causes the motion, dynamo action in the core 
transforms kinetic energy of the motion into magnetic energy of the field. Some of this 
energy emerges as our largely dipolar magnetic shield. But sometimes the transformation 
might go the other way, with the core running like a motor instead of a generator. This 
may accelerate core motions, weaken our shield, and perhaps drive our main axial dipole 
field to zero and beyond - to a reversed polarity. 

Paleomagnetic studies of crustal rocks magnetized in the geologic past reveal that polarity 
reversals have occurred many times during Earth’s history. Networked super-computer 
simulations of core field and flow, including effects of gravitational, pressure, rotational 
Coriolis, magnetic and viscous forces, suggest how this might happen in detail. And 


space-based measurements of the real, time-varying magnetic field help constrain 
estimates of the speed and direction of fluid iron flowing near the top of the core and 
enable tests of some hypotheses about such flow. 

Many core flow estimates omit narrow scale features, less than about 600 miles across. 
This is in part because observed changes in the main field are broad scale, so narrow scale 
flow is not easily resolved. Some estimates also assume the magnetic forces near the top 
of the core are negligible compared with other forces. Of course, doubts abound 
concerning both the scale of the flow and the strength of magnetic forces in a liquid metal 
buried beneath 1800 miles of rock. 

Now scientists at NASA’s Goddard Space Flight Center have developed and applied 
methods to test the hypotheses of narrow scale flow and of a dynamically weak magnetic 
field near the top of Earth’s core. Using two completely different methods, C. V. 
Voorhies has shown these hypotheses lead to specific theoretical forms for the 
“spectrum” of Earth’s main magnetic field and the spectrum of its rate of change. Much 
as solar physicists use a prism to separate sunlight into its spectrum, from long 
wavelength red to short wavelength blue light, geophysicists use a digital prism, spherical 
harmonic analysis, to separate the measured geomagnetic field into its spectrum, from 
long to short wavelength fields. They do this for the rate of change of the field as well. 

To test a hypothesis, Voorhies fits a theoretical spectrum to an observational spectrum. 
This yields an estimate of the radius of Earth’s core and its uncertainty. If the estimate 
agrees with the well-established core radius, determined by independent seismological 
techniques, then the hypothesis passes the test. If not, it fails. The weak field and narrow 
scale flow hypotheses each pass key tests, both separately and together. Recent analysis 
of data from the Danish Geomagnetic Research Satellite, however, suggests the core flow 
may have gotten a bit broader in scale, or better organized, over the past 20 years. 

For further information, including how compact eddies in a fluid conductor cause 
magnetic field changes over all scales, and how the magneto-geostrophic vorticity 
balance constrains kinetic-to-magnetic energy conversion atop the geodynamo, see 
“Narrow scale flow and a weak field by the top of Earth’s core: evidence from 0rsted, 
Magsat, and secular variation” by Coerte V. Voorhies, submitted for publication in the 
American Geophysical Union’s Journal of Geophysical Research - Solid Earth. 
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Abstract: To test two hypotheses against seismology, the 0rsted Initial Field Model is used to 
estimate the radius of Earth’s core by spectral methods. The model coefficients are used to 
compute the mean square magnetic flux density in spherical harmonics of degree n on the 
reference sphere of radius a = 6371.2 km, which is an observational spectrum R n . The theoretical 
spectrum tested, { R n c } = K(n+ll2)[(n(n+l)] A (d a) 2n+A , is obtained from the hypotheses of narrow 
scale flow and a dynamically weak magnetic field near the top of Earth’s core; it describes a low 
degree, core-source magnetic energy range. Core radius c and amplitude K are estimated by 
fitting log-theoretical to log-observational spectra at low degrees. Estimates of c from /? n of 
degrees 1 through N vary between 3441 and 3542 km as N increases from 4 to 12. None of these 
estimates differ significantly from the seismologic core radius of 3480 km. Significant 
differences do occur if N exceeds 12, which is consistent with appreciable non-core, crustal 
source fields at degrees 13 and above, or if other spectral forms are assumed. Similar results are 
obtained from the 1980 epoch Magsat model CM3. One way to deduce {Rf} uses an 
expectation spectrum for low degree secular variation (SV) induced by narrow scale flow near 
the top of Earth’s core, {F n c } = Cra(«+l/2)(«+l)(c/<2) 2n+4 The value of c obtained by fitting this 
form to the mean observational SV spectrum from model GSFC 9/80 is 3470 + 91 km, also in 
accord with seismologic estimates. This test of the narrow scale flow hypothesis is independent 
of the weak field hypothesis. The agreement between SV, Magsat, 0rsted and seismologic 
estimates of core radius means the hypotheses pass these tests. Additional tests are described. 

1. Introduction 

Among many physical hypotheses about Earth’s interior, let us develop and test two of 
geomagnetic and geodynamic interest. The first is the hypothesis of narrow scale fluid flow by 
the top of Earth’s core, which concerns core kinematics and geomagnetic secular variation (SV; 
see, e.g., Roberts & Scott [1965], Backus [1968]). “Narrow scale” here includes horizontal 
length scales less than about 10 3 km, such as the 4.2 km scale of Benton [1992], Many 
inversions of global geomagnetic change for core surface flow estimate only a broad-scale flow 
with features over 10 3 km across (see, e.g., Voorhies [1986, 1993, 1995]). 

The second hypothesis is that of a dynamically weak Lorentz force, or weak magnetic field, 
by the top of Earth’s core. This concerns core dynamics and the likelihood that the deep mantle 
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is too poor a conductor to allow much electric current across the core-mantle boundary (CMB; 
see, e.g., LeMouel [1984], Gubbins & Roberts [1987], Voorhies [1991]). “Dynamically weak” 
here includes force magnitudes much less than those of gravity, pressure and the Coriolis effect, 
albeit not necessarily as weak as the net force driving broad-scale accelerations of about 2xl0' 13 
m/s 2 [Voorhies, 1995]. This surficially weak field hypothesis is, however, compatible with a 
strong toroidal magnetic field deeper in the core (see, e.g., Backus [1986], Roberts & Gubbins 
[1987], Benton [1992]). A test might confirm that the potential field on Earth’s surface mainly 
represents a poloidal field from within the core, instead of from a strong toroidal core field 
threading CMB topography - hence a slight bending of otherwise toroidal field and poloidal 
coupling currents at a somewhat smaller jump in conductivity across the CMB. 

There are alternative hypotheses which have been tested; yet close fits to geomagnetic change 
obtained with broad-scale, non-geostrophic, core surface flows do not necessarily imply that 
narrow scale flow and/or weak field hypotheses are in error (see, e.g., Voorhies [1995]). 
Similarly, presence of some narrow scale flow does not imply absence of all broad-scale flow. 

Our hypotheses are embedded in a larger theory and a vast body of geomagnetic observations 
which ease their conversion into quantitatively testable forms. To this end, let us establish 
notation and a CMB layer model in sections 2 & 3; develop and test the narrow scale flow 
hypothesis in sections 4 & 5; develop and test the weak field hypothesis in sections 6 & 7; and 
offer a constraint on either the conductivity or the magnetic Reynolds number of Earth’s core. 

2. Notation and Background 

Let B(r, t) denote magnetic flux density at time t and position r in geocentric spherical polar 
coordinates (r, 0, (j>) caused by electric current of density J and magnetization M within the earth. 
Above Earth’s surface, solenoidal B equals the negative gradient of the zero-mean scalar internal 
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magnetic potential V. This potential satisfies Laplace’s equation and has a Schmidt-normalized 
spherical harmonic expansion with Gauss coefficients of degree n and order m, denoted [g n m (t), 
h n m (t)] on a reference sphere of radius a = 6371.2 km. Coefficients through finite degree TVp can 
be estimated by spherical harmonic analysis of the measured field, which also enables separation 
of the internal-source field considered here from external-source fields (see, e.g., Langel [1987]). 

As is also well-known, the mean square field represented by harmonics of degree n, averaged 
over a sphere of radius r containing the sources, is given by 
n 

Rn(r, t) = (n + l)(o/r) 2n+4 2 [g n m (t)f + [h n m (t)] 2 (1) 

m = 0 

(see, e.g., Lowes [1966], Lucke [1957], Mauersberger [1956]; Meyer [1985]). The /? n form a 
discrete geomagnetic spectrum with units of T 2 . The integral of magnetic energy density over the 
sphere has a related spectrum, 2nr 2 RJ\i 0 in vacuum permeability |i 0 , with units of spatial power 
or force (J/m). A spherical harmonic expansion is equivalent to a centered, 2n-pole moment 
expansion, so each R n represents a multipole power. 

The mean square rate of change of the magnetic field represented by harmonics of degree n, 
averaged over the same sphere, is given by the S V spectrum 
n 

F a (r,t) = (n + l)(a/r) lD+4 I [0 t g n m ) 2 + OtCV] • (2) 

m=0 

Note that F n is not 9 t R D ; moreover, the time average of magnetic spectrum (1), denoted <R a >, is 
not generally the spectrum of the time averaged field. Similarly, when SV itself varies in time, 
the time average of SV spectrum (2), denoted <F n >, is not the spectrum of net magnetic change. 

Spectra computed from coefficients determined by weighted least squares fits to geomagnetic 
measurements, with no assumptions about a core field, are here called “observational spectra”. 
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The accuracy of an observational spectrum depends on geomagnetic field modeling procedure, 
such as the choice of Afc, as well as data selection, distribution and accuracy. Spectra obtained 
from physical hypotheses about sources are here called “theoretical spectra” (see, e.g., McLeod 
[1996]; Voorhies [1998]; Voorhies, Sabaka & Purucker [2002]). Hypotheses about averages 
over physical processes may yield a theoretical spectrum that is also a mean, or “expectation 
spectrum”, denoted {/?„} for the main field or {F„} for the SV. 

Observational spectra R„ determined by analyses of satellite magnetic data have been 
interpreted in terms of a mainly core-source field for n < 14 and a dominantly crustal -source field 
for n > 14 (see, e.g., Langel & Estes [1982], Cain et al. [1989], Voorhies, Sabaka & Purucker 
[2002]). A review of publications bearing on a core-source interpretation for portions of R n and 
F n is omitted for brevity (see, e.g., Booker [1969]; Verosub & Cox [1971]; Lowes [1974]; 
McLeod & Coleman [1980]; Hide & Malin [1981]; Langel & Estes [1982]; Benton et al. [1982]; 
Shure, Parker & Backus [1982]; Voorhies & Benton [1982]; Gubbins [1983]; Stevenson [1983]; 
Voorhies [1984]; Gubbins & Bloxham [1985]; McLeod [1985, 1996]; Meyer [1985]; Benton & 
Alldredge [1987]; Benton & Voorhies [1987]; Backus [1988]; Constable & Parker [1988]; Cain 
et al. [1989]; Hulot, LeMouel & Wahr [1992]; Harrison [1994]; Hulot & LeMouel [1994]; 
Voorhies & Conrad [1996]; Walker & Backus [1997]; Voorhies et al. [2002]; De Santis, 
Barraclough & Tozzi [2003]). 

Clearly, an observational spectrum may be used to help test a theoretical spectrum, hence its 
underlying physical hypotheses. There is, however, a third class of spectra, here called 
“constrained spectra”, computed from coefficients constrained by one or more assumptions about 
a core field (see, e.g., Shure, Parker & Backus [1982]; Gubbins [1983]; Backus [1988]). The 
independence, hence utility, of a constrained spectrum can be compromised by constraints that 
force it to either agree or disagree with a theoretical spectrum. In special circumstances. 
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however, the statistical significance of increased misfit to measured data caused by a constraint 
could provide a test of physical hypotheses underlying the constraint. 

Theoretical spectra from sources in a roughly spherical core of radius c are here denoted /? n c . 
For example, Gubbins’ [1975] expression for the minimum value of the finite Ohmic dissipation 
in the core implies that, for degrees in a magnetic dissipation range defined by n > Afo, 

R n c (a- n > N d ) < Kg m 2 5 (c/a) 2n , (3a) 

where constant K G > 0 and 5 > 0. Though finite, N D may be much larger than 12. So there may 
be a low degree magnetic energy range, defined by n < Ns < N D , and perhaps an intermediate, if 
not inertial, sub-range defined by Ns < n < This is the case for theoretical spectra of 
Stevenson [1983] and McLeod [1985, 1996], as shown in Appendix A. A scale analysis 
illustrates the physical plausibility of a generalized Stevenson - McLeod relation, 

{R u c (r> c-n<N E )} » K (n + Viftnin+l)]' 1 (c/r) 2n+4 , (3b) 

as an expectation spectrum for low degrees, as shown in Appendix B. 

Both the scale analysis and the empirical approach in section 6 yield spectrum (3b) as an 
expected consequence of narrow scale flow and a weak field by the top of the core. Neither 
specify expected spectral variance, denoted {(R n c - {A n c }) 2 }; therefore, tests of spectrum (3b) are 
here limited to comparisons between magnetic estimates of c and independent estimates. 
Seismologic estimates of c differ by but a few km and are here denoted c s = 3480 km (see, e.g., 
Dziewonski & Anderson [1981]; Kennett, Engdahl & Bulland [1995]). 

3. CMB Layer and Narrow Scale Flow 

The transition from fluid conducting core to rigid resistive mantle has long been modeled 
simply by a sharp, impenetrable interface: a fixed jump in material properties across the CMB. 
Finite Lorentz and viscous forces are ensured by continuity of both magnetic and hydrodynamic 
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stress tensors across this interface. Then B is also continuous and the relative velocity u of a 
Newtonian fluid that wets the mantle satisfies both kinematic and no-slip CMB conditions. 

Fluid motion at depth indicates a boundary layer between the interface and a main stream. 
Analyses of quasi-steady magnetic, mass, and momentum transport equations indicate a thin, 
weak, boundary layer that neither generates much electrical current nor absorbs appreciable 
normal fluid flow (see, e.g., Ball, Kahle & Vestine [1969], Hide & Stewartson [1972], Benton 
[1981], Gubbin's & Roberts [1987], Benton [1992]). This depends on a kinematic shear viscosity 
v that is very small and a magnetic diffusivity r\ that is not too large, conditions met with v ~ 
3xl0' 7 m 2 /s and r\ = (poo)" 1 = 1.6 m 2 /s for electric conductivity a = 5x10 s S/m [Poirier, 1988; 
Lumb & Aldridge, 1991; Voorhies, 1999; Dobson et al., 2000], A thin viscous sub-layer is also 
suggested by equating typical magnitudes of tangential viscous forces with Coriolis or Lorentz 
forces and solving for Ekman or Hartmann scale depths, 5g = 8 cm or S H = 15 cm, respectively. 
The induction equation for uniform c, our non-conservative magnetic transport equation 
d,B = Vx(u x B) + qVxVxB , (4) 

shows core-source SV can be induced by mainly lateral motion at the base of a viscous sub-layer. 
A feeble boundary current implies a small jump in B across the sub-layer, negligibly small for the 
normal component. Though purely diffusive at the interface itself, the radial component of (4) 
has thus been used to analyze SV in terms of broad-scale fluid flow, and occasionally flux 
diffusion, at the top of a spherical main stream. Constrained inversions of geomagnetic secular 
change indicate a typical flow speed U of about 7.5 km/yr (see, e.g., Voorhies [1995]). 

With the foregoing values for U, 5e, v and c s , the main stream Reynolds number Uc s /v of 
about 3xl0 9 is so much greater than the boundary number U8s/v of about 60 as to indicate some 
small scale motions near the top of the core, perhaps in a thicker, unsteady, second boundary 
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layer featuring eddy mixing and enhanced diffusion - if not entrainment of the overlying viscous 
sub-layer. Instead of a typical width for fronts between broad regions of more uniform flow, the 
4.2 km lateral length scale for fluid velocity obtained by Benton [1992], denoted k, may describe 
a seething mass of short-lived, rotationally polarized hydromagnetic eddies. By equation (4), 
such narrow scale eddies could individually induce narrow scale field variations, yet could also 
contribute collectively to observable, broad-scale SV. 

A quantitative model of this contribution is needed to test the narrow scale flow hypothesis. 
A model consisting of pseudo-random walks of magnetic field line foot points with a single eddy 
diffusivity Ulo ~ 1 m 2 /s seems too much like molecular magnetic diffusion to test; moreover, a 
linear diffusion term with scale-invariant diffusivity cannot describe how narrow scale flow 
induces broad-scale SV by mode-mixing. Deterministic inversions of (4) do not resolve narrow 
eddies and might misattribute eddy mixing of magnetic modes to broad-scale flow or flux 
diffusion. Deterministic forward models, notably numerical dynamo models that solve a system 
of magnetic, mass, momentum and energy transport equations closed by an equation of state (see, 
e.g., Glatzmaier & Roberts [1995a, b]), can better resolve compact eddies and effects of mode 
mixing. 

For example, a dynamo simulation used to investigate the frozen-flux core approximation 
included degrees as high as 239 [Roberts & Glatzmaier , 2000]; yet observational field models 
similarly used included degrees of at most 13 and often 10 or less [Hide & Malin, 1981; Voorhies 
& Benton, 1982; Voorhies, 1984; Benton & Voorhies, 1987]. With l = [ 87 tc s 2 /n(n+l)] 1/2 , we find 
Z sim > 73 km = 17Zo and Z obs > 1300 km = 18Z sim ; therefore, finer resolution is needed to simulate 
eddies of scale l 0 . Additional assumptions about material properties, turbulent diffusivities, 
boundary conditions, and initial conditions further complicate hypothesis testing via numerical 
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simulation (see, e.g., Glatzmaier [2002]). Both forward and inverse deterministic models of SV 
induced by narrow scale flow are bypassed here via a statistical kinematic model of SV induced 
by compact eddies. 

4. An Equivalent SV Spectrum from Narrow Scale Flow 

On the narrow scale flow hypothesis, magnetic changes at different locations at the top of the 
core are induced by different eddies transporting different local fields in different directions at 
different speeds. Such changes may well appear uncorrelated on scales broader than the eddies. 
Mathematically, the magnetic change induced by each narrow scale eddy can be approximated in 
the far field, notably well outside the core, by an equivalent source of change at the base of a 
viscous sub-layer. Following McLeod [1996], quasi-static lateral magnetic transport during 
differential interval At causes differential exterior magnetic change AB equivalent to differential 
dipole moment changes Ad; scattered atop the main stream at radius c (i = 1, 2, 3, . . . , I). 

To see this, recall that a single magnetic flux vector Bo dA at fixed position xo on the surface 
A of the source region acts as the point source of a dipole field with moment proportional to 
Bo dA. The magnetostatic field at position x due to this equivalent source is well-known (see, 
e.g., Jackson [1975, equation 5.64]). Infinitesimal quasi-static lateral displacement Ax of this 
single magnetic vector, with no change in orientation and magnitude, would cause a net change 
in the exterior field equivalent to a differential quadrupole moment at Xo + Ax/2. More generally, 
however, there is a magnetic flux vector at each position on the source surface; lateral transport 
replaces the vector at x 0 with an adjacent vector of slightly different orientation and magnitude; 
and the change in the exterior field is equivalent to that of a differential dipole moment Ado at xo. 
Given many equivalent source changes Adj at x, on A (i = 1, 2, 3, ...,/), the total change in the 
exterior field AB at |x| > c follows by superposition. In the continuum, a differential change in 
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surface magnetic moment density at x’ on c replaces discrete Ad; as the equivalent source for 
exterior secular change AB(x,f), which follows by integration over x’. 

Elements of the dyad formed by two differential dipole moment changes at two well 
separated points, [Adi][Ad 2 ] T , may be either positive or negative. The average over a 
kinematically unbiased ensemble of such dyads gives zero cross-correlations, but non-zero auto- 
correlations. The resulting expectation spectrum for broad-scale SV is equivalent to that from 
laterally uncorrelated, randomly varying dipole moments on the source shell c, which is given by 

{F n c (r>c)} = Cn(n + Vi)(n + l)(cVr) 2n+4 . (5) 

Equation (5) differs slightly from McLeod [1996, equation (11)] because changes in horizontal as 
well as radial components of an equivalent source can contribute to {F n c }: for example, consider 
rotation about the vertical of a horizontal equivalent source at the equator. Positive amplitude C 
is proportional to {/(Adi/Af) 2 }. Physically, C tends to increase with transport speed, field 
gradients and intensity; yet the mathematical derivation of spectrum (5) in Appendix C does not 
use equation (4). It uses steps analogous to those in Voorhies [1998, equations (6a)-(20b)]. 

Laterally uncorrelated dipole changes offer a rough model of SV rich in narrow features; yet 
the sum of attenuated cubic S V spectrum (5) over n converges on all spheres of radius r>c>c. 
The cubic polynomial, which modulates the exponential attenuation (c'/r) 2n+4 of a potential field, 
increases with n faster than does the linear polynomial obtained from SV sources equivalent to 
uncorrelated Dirac delta-functions in 2? r . The latter are sources for the famous “white noise” 
spectrum, so SV spectrum (5) is said to be “blue” or “hard” - to borrow descriptors of spectra 
rich in short wavelength components from optical or X-ray spectroscopy, respectively. Of 
course, no single Ad represents perfectly the magnetic change induced by an extended eddy, so 
spectrum (5) will not represent SV on scales as narrow as the eddies themselves - hence at high 
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degrees (e.g., n > 4,200 for Z 0 above). Moreover, even transient cross-correlations may result in 
deviations from spectrum (5) at lower degrees. Indeed, a softer or less blue SV spectrum, such as 
an attenuated quadratic, may indicate contributions from partially resolved eddies. 

Each set of Adi implies one F n c , so derivation of theoretical SV spectral variance {[F n c - 
{F n c }] 2 } would require additional assumptions about the ensemble of equivalent SV sources - 
hence additional physical hypotheses about the eddy transport they represent. Extra hypotheses 
tend to complicate tests of the narrow scale flow hypothesis. Low degree physical deviations F n c 
- {F n c } caused by core dynamic processes can, however, be represented via cross-correlated 
equivalent SV sources. And cross-correlation of even a small fraction I' m of the Adj can cause 
deviations with magnitudes similar to {F n c }. We thus anticipate a spectral variance of magnitude 
similar to, or perhaps a few times larger than, { F n c } 2 itself. Other aspects of SV spectral 
covariance are noted in Appendix D. Spectrum (5) is arguably better tested against a time 
averaged observational spectrum, not only to help average out effects of transient cross- 
correlations, but to identify persistent deviations from it. 

5. Tests of the SV Spectrum 

Theoretical SV spectrum (5) has two parameters, amplitude C and source radius c, to be 
estimated by fitting observational SV spectra F n (a,t). For a thin sub-layer, an estimate of c 
amounts to an estimate of c, so the significance of its difference from seismologic core radius c s 
provides a test of equation (5), hence the underlying narrow scale flow hypothesis. Curiously, 
the misfit between theoretical and observational spectra does not provide a sensitive test of the 
hypothesis. This is because the theoretical spectral covariance needed to fully establish the 
statistical significance of such misfit is not specified (see Appendix E). We therefore emphasize 
comparison of estimated SV source radius c with independently determined c s . 
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To ease this estimation, and anticipating fluctuations about (5) amounting to a factor of about 
er x , we minimize the sum of squared residuals to observational ln(F n ) for degrees n^ n to n max . 
The sum of squared residuals per degree of freedom for d = n max - nmin + 1 data fitted by p 
parameters is just 

ttmax 

q 2 = (d-pf Z [ln(F D )- ln{F a c }) 2 . (6) 

^ = ^min 

Of course, only if the residuals were approximately log-normally distributed might one hope that 
the estimates approximate maximum likelihood estimates. The estimation requires computation 
of expected parameter covariance, with the square root of the diagonal variances indicating 
expected parameter uncertainties. Multiplication of these values by q yields scaled parameter 
uncertainty estimates. Typically q < 1, so the scaled uncertainties are less, and the sensitivity of 
the test greater, than expected. 

5.1 Initial Test. The observational SV spectrum fitted first was computed from field model 
of Langel, Estes & Mead [1982], This model GSFC 9/80 closely fits 15,206 Magsat 
observations, 71,000 POGO observations, measurements from 148 observatories, 300 filtered 
marine data, and 600 measurements from select repeat stations. It includes main field, first, 
second, and third time derivative coefficients through degrees 13, 13, 6, and 4, respectively, and 
observatory biases to account for local lithospheric magnetic anomalies. Unlike some more 
recent field models, it imposes no smoothness or other constraints upon a core-source field. 

Table 1 lists epoch t of F a (a,t), the range of degrees n fitted (n^n-nmax), estimated core radius 
c" and its scaled uncertainty for three ranges. The bottom line, labeled “Avg.”, gives results from 
fitting theoretical spectrum (5) to <F n > integrated over the interval 1960-1980 spanned by the 
model. The fit to all degrees 1-13 of the time-averaged spectrum yields c = 3470 + 91 km. Table 
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1 also shows results from fitting intermediate degrees 3-11 and, due to long-standing concerns 
about SV coefficients above degree 10 [Voorhies, 1984], from degrees 3-10. 

The agreement between geomagnetic SV spectral estimates of core radius in Table 1 and the 
3480 km seismologic value is excellent. Indeed, only one value of 18 differs from c s by more 
than twice its scaled uncertainty. The tabulated results are summarized as c = 3.5 + 0.1 Mm = c s . 
The hypothesis of narrow scale flow by the top of Earth’s core, as represented by the broad-scale 
SV spectrum (5) expected from an ensemble of compact eddies inducing laterally uncorrelated 
S V, thus passes our initial test. 

For ft min = 1, however, closer study reveals a systematic increase in estimates of c from 
GSFC 9/80 with ft max . This is shown in Table 2a, which lists degree range, q from equation (6), 
estimated core radius with scaled uncertainty estimate in km, and error relative to c s in km. The 
latter are judged significant for n max <11 and so indicate a failure of theoretical spectrum (5); 
moreover, if the larger misfits found for ftmax > 10 are due to poorly determined SV coefficients, 
then the smaller errors might be fortuitous. To check this, we set ftmi n = 3; as shown in Table 2b, 
so doing tightens the fit and eliminates errors in c in excess of twice the scaled uncertainties. 
Table 2b also shows increased misfit from degrees above 10. The tabulated results imply 
theoretical spectrum (5) adequately describes this observational SV spectrum except for degrees 
one and two. Downward continuation shows these exceptional terms contribute little to SV at 
the CMB. They may be reconciled via degree-dependent process variances from an alternative, 
non-log-normal, distribution. The exceptions can be understood in terms of fast decline of a 
strong dipole and rapid quadrupole rebound as defined below. 

Voorhies & Conrad [1996] found R\ greater, and R 2 less, than expected based on fits of 
spectra (Ala) or (Alb) to observational spectra at epoch 1980 - albeit within the ranges we 
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expected 80% of the time. We also found R\ decreasing and Ri increasing; moreover, for all 
orders m, we found (5tgi m )/gi m < 0 and (9^2™ )lg 2 n > 0. The chances of such perfect (anti-) 
correlations were put at 1/8 for the dipole, 1/32 for the quadrupole, and 1/256 jointly. This 
otherwise remarkable coincidence was viewed merely as an efficient relaxation of the field 
towards expectation values. The large value of F 2 /R 2 and the perfect correlation, and the large 
value of dxR^Ri itself, are called “rapid quadrupole rebound”, as distinct from the quadrupole 
diminution noted by Stevenson [1983]. 

5.2 A Softer SV Spectrum? To further investigate SV spectra, we use the attenuated 
quadratic SV spectrum expected from uncorrelated, randomly varying dipoles scattered 
throughout the interior of a ball of radius c, which is given by 

{F u c (a>e)}* = C*n(n + l)(c/a) 2n+4 . (7) 

This softer SV spectrum is thought to be a proxy for contributions from some partially resolved 
eddies. It is not thought to describe SV from small scale eddies scattered throughout the core, or 
even a layer thicker than about 90 km, because equation (4) fixes the origin of core-source SV at 
the top of the core. Surprisingly, spectrum (7) may also be a proxy for effects of laterally 
heterogeneous mantle conductivity. Laterally homogeneous conductivity tends to harden a core- 
source SV spectrum because physical attenuation decreases with harmonic degree [ McDonald , 
1957]; however, mode coupling by lateral heterogeneity in deep mantle conductivity may in 
effect scatter some intense, narrow scale, core-source SV into broader scales, thereby softening a 
core-source SV spectrum before it emerges through Earth’s surface. 

Tables 3a and 3b are analogous to Tables 2a and 2b, but show results of fitting proxy 
spectrum (7) to the 20 year average <F n > from model GSFC 9/80. Table 3a, with = 1, shows 
small and insignificant errors in c for <11. Table 3b, with = 3, shows large positive 
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errors in c that exceed twice the scaled uncertainty for R mi > 8. Evidently, spectrum (7) 
accommodates rapid dipole decline and quadrupole rebound, but is too soft for higher degrees. 
Further tests against other observational spectra might better distinguish between spectra (5), (7), 
and promising intermediate modulation factors such as [n(n+l)] 5/4 , particularly if the maximum 
degree of reliable F n (t) can be established. 

5.3 More Tests. Model CM3 [Sabaka, Olsen & Langel, 2002], fitted to Magsat, POGO and 
observatory data from 1960 to 1985, features a more comprehensive representation of external 
source fields than earlier models, an internal static field through degree 65, and a temporal spline 
parameterization of SV. The SV model was, however, constrained to reduce the amplitude of 
narrow scale SV that is poorly determined by geographically sparse data before POGO and after 
Magsat satellite surveys. Two constraints were used. The first conflicts with spectrum (5) by 
forcing the mean square value of the surface Laplacian of 3B r /3t, averaged over the sphere of 
radius c s and time, to be small - as if the SV spectral modulation factor were n 4 or less instead of 
n +3 . This first constraint was not as strongly imposed as the second, which forces the mean 
square second time derivative of B r , also averaged over the sphere of radius over c s and time, to 
be small. This reduces temporal variability in F n from model CM3, notably at high degrees; 
indeed, it helps makes the F n steady to within 18% for n > 5, and to within 5% for n > 10. 

Comparison of mean SV spectra from models CM3 and GSFC 9/80, both averaged from 
1960-1980, shows the two sets of <F n > agree to within 25% for n < 11, with CM3 giving smaller 
values - typically 12% smaller. For degrees 11, 12, and 13, however, values from GSFC 9/80 
exceed those from CM3 by factors of 2.6, 1.8, and 34.5, respectively. The latter exceeds the 
factor of e anticipated from process variance. Evidently, <Fi 3 > is not reliably determined and so 
is not considered further. 
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Tables 4a and 4b show results of fitting theoretical S Y spectrum (5) to the mean SV spectrum 
from CM3, averaged from 1960 to 1980. Low values for c in Table 4a might again suggest 
spectrum (5) is too hard; yet the increase in estimated c with n max is evident. Table 4b confirms 
this to be largely due to <Fj> and <F 2 >. The fits listed in Table 4b are quite tight; indeed, for the 
minimum q 2 of 0.0321, deviations from {F n } are typically a factor of (1 .2)— 1 instead of er. 

Tables 5a and 5b show fits of proxy SV spectrum (7) to CM3. They confirm that the 
attenuated quadratic form does fairly well considering degrees 1 through n max < 13, but yields 
errors in c that are judged significant for higher degrees 3 though n max > 7. Downward 
continuation shows the higher degrees contribute far more to SV by the CMB than do F\ and F 2 . 
Although extended study suggests proxy spectrum (7) might be more suitable during another, 
shorter epoch (see Appendix F), it is bypassed for now in favor of expectation SV spectrum (5). 

6. Core Magnetic Spectrum from Narrow Scale Flow and a Weak Field 

Though (F a /R n ) was studied empirically as a kind of summary dispersion relation for the 
core field, M. G. McLeod (1985, pers. comm.) pointed out other reasons to consider spectral ratio 
Rn/F n c . As it turns out, some hypotheses indirectly constrain this ratio in ways that allow a 
theoretical core-source spectrum R n c to be obtained from a theoretical SV spectrum F n c . One 
such hypothesis specifies the form of the temporal power spectrum [McLeod, 1996]. Another is 
the hypothesis of “constant aspect ratio”, which asserts a single direct proportionality between 
horizontal wave-numbers and effective radial wave-numbers in the non-potential portion of the 
poloidal field near a spherical CMB. Such radial wave-numbers arise from analysis of the radial 
component of induction equation (4) at the top of a viscous sub-layer. A third also considers the 
radial component of (4), but at the top of a free-stream where narrow scale eddies mix and remix 
modes of a dynamically weak magnetic field. Yet another reorders arguments in Appendix B.3 

to obtain a relation proportional to (B17), hence (B21b), from relations (Bllb), (B18) and (4). 
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All these hypotheses lead from attenuated cubic SV spectrum (5) to attenuated 1 In core field 
spectra like (3b) and all seem compatible with the weak field hypothesis. For brevity, it is here 
argued on dimensional grounds that the available diffusivities, the weak field hypothesis, and SV 
spectrum (5) lead to this form for {R n c } . 

Recall our narrow eddies at the base of a viscous sub-layer. Each eddy (i = 1, 2, 3, ... 7) has a 
characteristic lateral speed C7(i), a characteristic lateral length scale L(i) « c, hence a lateral eddy 
diffusivity U(i)L(i). The ratios of lateral eddy diffusivities to molecular magnetic diffusivity r) 
define eddy magnetic Reynolds numbers U(i)L(i)/r\ = A(i). By continuity, the normal component 
of eddy-induced and other SV crosses a thin, weak viscous sub-layer largely unaltered and 
specifies the SV signal emerging from a spherical core into a source-free exterior. At the top of 
the sub-layer, however, u vanishes and even the eddy-induced portion of the signal is transmitted 
via molecular diffusion. Owing to its importance in such radial magnetic transport, the 
molecular diffusivity may be considered a radial diffusivity. The A(i) may then be considered 
indices of diffusive anisotropy - ratios of lateral eddy to radial molecular diffusivities. 

Neither r| nor eddy diffusivities are directly measured; however, observation and analysis can 
reveal empirical diffusivities. In particular, for a sufficiently time-varying field, time-averaged 
observational spectra </?„> and <F n > together define regular empirical time constants 

r n ^ [</? n >/<F n >] 1/2 . (8) 

When combined with horizontal wave-numbers, defined via the surface Laplacian operator to be 
O = [«(n+l)/c 2 ] 1/2 , these T n further define lateral empirical diffusivities 

D a S (k Q h y 2 T n ' = [c 2 /n(n+ 1 )] [<F a >/<R n >] 1/2 . (9a) 

- The theoretical counterpart to definition (9a) is written 

Cn = (kn)V = [c 2 / u(n+ 1 )] [ { F n c } / { R n c } ] 1/2 . (9b) 
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For a core-source field governed by induction equation (4), the only physical diffusivities that 
can depend on are r| and the eddy diffusivities, though the latter may depend on other quantities. 
On dimensional grounds, any dependence of on degree n should, and arguably must, be 
determined from that of these diffusivities. 

The magnetic field and Ohmic heating of interest are far too weak to cause either appreciable 
anisotropy or heterogeneity in core electric conductivity; therefore, r) is effectively independent 
of the field, its SV and the harmonic degrees thereof. Eddy speeds and length scales, hence eddy 
diffusivities, do not depend directly on SV. They depend only on the fluid velocity, which can 
depend on the field via the Lorentz force, JxB, in the usual momentum transport equation. 

If the Lorentz force is weak compared with other forces near the core surface, as for 
tangentially geostrophic flow [LeMouel, 1984], then the fluid velocity and the eddy diffusivities 
depend but weakly on the magnetic field. If this dependence is negligible, then the eddy 
diffusivities must be effectively independent of the magnetic field, hence the harmonic degrees of 
the field. On this weak field hypothesis, neither molecular nor eddy diffusivities near the top of 
the core depend on the harmonic degrees n of either the core field or its SV. Yet these are the 
physical diffusivities from which the dependence of C, D on degree n must be determined. With no 
basis on which to construct a degree-dependent theory of £ n , we deduce that the in (9b) must 
reduce to a single constant, independent of n, 

Cn = [c 2 /n(n+l)][{F n c }/{i? n c }] 1/2 = C, (10a) 

or 

{F n C }/{*n C } = CWn+1)] 2 = * 2 . (10b) 

By (10b), the core field spectrum is expected to be much softer than that of core-source SV. 
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To check (10b), Voorhies & Conrad [1996] fitted the function fioln[n(n+l)] - Undo to 
observational values of ln[F n /R n ] from degrees 3-12 of model GSFC 9/80 at epochs 1960, 1970, 
and 1980. The three resulting values for (3 0 average to 1.957 ± 0.156. This agrees with the 
expectation value (3 0 = 2 in (10b). The check is independent of the radius of Earth’s core. The 
three values for oto are within a factor of 1.94 of 2,640 years (not the “26,400” years misprinted in 
Voorhies & Conrad [1996]). The implications for £, R m , and a are discussed in Appendix G. 

Substitution of expectation SV spectrum (5), from the narrow scale flow hypothesis, into 
equation (10b), from the weak field hypothesis, yields our expectation core field spectrum 


{/? n c } = C£V(n + Vi)[n{n+ l)]" 1 (c/ r) 2n+4 

(Ha) 

= K(n + Vi)[n(n+ 1 )]"’ (cVr) 2n+4 . 

(lib) 


Though obtained in different ways, quantitative distinctions between spectrum (1 lb) and earlier 
forms (A.la), (A.lb), and (A.lc) are largely confined to degrees 1 and 2. 

7. Tests of the Core Field Spectrum 

Expectation spectrum (lib) has two parameters, amplitude K and source radius c, to be 
estimated by fitting observational spectra R n (a,t). An estimate of c again amounts to an estimate 
of c, so the significance of its difference from c s provides a test of (lib), hence the underlying 
hypotheses of narrow scale flow and a dynamically weak field by the top of Earth’s core. The 
hypotheses do not specify process variance { [i? n c - {i? n c }] 2 } , so the misfit between theoretical and 
observational spectra does not provide a sensitive test of the hypotheses. This process variance 
might be large; therefore, we again emphasize comparison of magneto-spectral estimates of c 
with the seismologic value c s . 
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To ease this estimation, and anticipating fluctuations about (lib) amounting to a factor of 
about e ±1 , we minimize the sum of squared residuals to observational ln(R n ) for degrees n^n to 
ttmax- The sum of squared residuals per degree of freedom for n ^ - n m i„ + 1 data fitted by 2 
parameters is just 

^inax 

S 2 = (n max - rimin - 1)’ 1 2 [InfRJ- ln{R n c }] 2 . (12) 

^ = ^min 

The estimation requires computation of expected parameter covariance; the square root of the 
variances give expected parameter uncertainties. Multiplication of these values by s yields scaled 
parameter uncertainty estimates. Typically, s < 1 and the test is more sensitive than expected. 

Voorhies et al. [2002] describe a test of spectrum (lib) at Magsat epoch 1980. Here 
spectrum (lib) is tested against the independent observational spectrum from the 0rsted Initial 
Field Model (OEFM) [Olsen, et al., 2000]. The epoch 2000 OIFM features a weighted least 
squares fit of main field coefficients through degree 19, and external field coefficients of degrees 
1 and 2, to 13,859 select data acquired by the Danish Geomagnetic Research Satellite 0rsted. 
Table 6 lists the range of degrees fitted (n mn = 1 though rc max ), s 2 , the estimate of c with scaled 
uncertainty, and the error of the estimate relative to c s = 3480 km. Selection of the minimum s 2 
solution at = 12 fixes a third parameter and yields c = 3542 + 61 km as the core radius 
estimated from 0rsted. No significant errors in c are found for degree ranges 1 to n max < 13, so 
(1 lb) passes this test. Estimates from these degree ranges all agree and average to 3489 + 39 km. 
The increased misfit and significant errors introduced with degrees above 12 are attributed to 
non -core, likely crustal, source fields. 

To refine the test by Voorhies et al. [2002], we use R n from model CM3 of Sabaka, Olsen 
and Langel [2002] at Magsat epoch 1980. Table 7 lists the range of degrees fitted, s 2 , the 
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estimate of c with scaled uncertainty, and the error of estimate relative to c s . Again, no 
significant errors in estimates of c are found for degree ranges 1 to n max <13. Radii from this 
degree range all agree and average to 3495 ± 28 km. This average more heavily weights lower 
degree /? n , unlike the process variance for the R n c distribution of Voorhies & Conrad [1996]. 

Table 8 is analogous to Table 7, but with = 3. It confirms that much of the scatter, as 
measured by the larger values for s 2 in Table 7, comes from strong Ri and weak Rj. Table 8 also 
shows the effect on estimates of c associated with R% being lower, and R 9 higher, than expected 
based on the fit of (lib). The larger residuals at degrees 1, 2, 8, and 9 may help provide some 
indication of geomagnetic variability about (lib), hence process variance. 

A reviewer asks why analysis of magnetic spectra is any better at determining core radius 
than the frozen-flux method of Hide. I did not claim it is; however, estimates obtained using 
spectrum (lib) and the main field models are typically more accurate and precise than those 
obtained using the frozen-flux approximation [Hide & Malin, 1981; Voorhies & Benton, 1982; 
Voorhies, 1984], In particular, some 44 frozen-flux core locations obtained from a few field 
models at various truncation levels average to 3506.2 + 300.9 [Voorhies, 1984, equation (3.21)]. 
The first 9 estimates in Table 6 average to 3489 + 39 km; the first 9 in Table 7 average to 3495 + 
28 km. Frozen-flux methods rely heavily upon uncertain phase information in harmonic orders m 
and upon uncertain secular change information from either SV models or main field models at 
different epochs. Analysis of F n also relies on SV. In contrast, the main field spectral method 
relies on comparatively well-determined R a alone; moreover, it does not require the frozen-flux 
approximation - either in section 6, in the scale analysis of Appendix B.3 and B.4 that requires 
finite o to obtain spectrum (3b), or in the analysis of flux diffusion that returns a similar 
spectmm on the constant aspect ratio hypothesis. Evidently, these advantages reduce scatter. 
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8. Summary and Conclusions 

A theoretical form for the low degree, core-source geomagnetic spectrum {R n c } is deduced 
from the hypotheses of narrow scale flow and a dynamically weak magnetic field by the top of 
Earth’s core. This form (lib) differs but slightly from those advanced by Stevenson [1983] and 
McLeod [1985, 1996]. To test these hypotheses, this theoretical spectrum is fitted to two 
observational spectra, one determined by analysis of data from the 0rsted satellite and the other 
by analysis of independent data from Magsat and POGO satellites and surface observatories. 
These fits yield estimates of the radius of Earth’s core. These estimates differ insignificantly 
from the seismologically established radius, as judged by the scaled uncertainty estimates, for 
degree ranges 1 to n max <13. Because these errors of estimate are not significant, the hypotheses 
taken together pass these tests. 

Significant errors introduced with higher degree multipole powers are attributed to non-core, 
likely crustal source fields. Very small errors, less than 10 km, are found by excluding the 
strong, rapidly declining dipole power and the weak, rapidly rebounding, quadrupole power. 

To deduce the expectation spectrum (lib), we can use the theoretical low degree SV 
spectrum (5), which follows from the narrow scale flow hypothesis alone. There are other ways 
to obtain both spectra (lib) and (5), as illustrated in Appendix B; however, the geophysical 
foundation of all such methods known to date appears to rest on the narrow scale flow and weak 
field hypotheses - provided the temporal power spectrum of McLeod [1996] can be inferred from 
the approximate spatial spectra. Initial tests of SV spectrum (5) gave estimates of the radius of 
Earth’s core that differ insignificantly from the seismologic value. This shows that the narrow 
scale flow hypothesis can pass a test against seismology that is independent of the surficially 
weak field hypothesis. Attempts to test (5) with the same rigor as (lib), however, reveal errors 
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arising from rapid dipole decline and quadrupole rebound. Though Fi and F 2 contribute little to 
SV at the CMB, such exceptional behavior of may cast doubt on (5); yet it may instead indicate a 
need to use degree-dependent process variances and a different distribution for residuals about 
SV spectrum (5). 

Clearly, physical hypotheses are not proved true simply by passing a few tests, nor does 
falsification preclude their use in successive approximation. And there might be different 
physical hypotheses that yield theoretical spectra like those presented here. The hypotheses of 
narrow scale flow by the top of the core, as described by expectation SV spectrum (5), and of a 
dynamically weak core surface field, described via magnetic spectrum (lib), have nonetheless 
demonstrated considerable merit. By enabling magneto-locations of Earth’s core, these spectra 
demonstrate greater utility and accuracy than far softer spectra. Evidently, far softer spectra do 
not provide reliable prior information on either the broad scale field or its secular variation. It is 
hoped these geophysical hypotheses and geomagnetic spectral forms will be further tested, and 
provide some basis for comparison with alternatives that may emerge, in the future. 

Appendix A: Magnetic Spectral Ranges 

To compare a core-source field spectrum with low degree energy, intermediate, and high 
degree dissipation ranges with theory, recall the magnetic energy spectrum M(k) as a function of 
Cartesian-Fourier wave-number k (see, e.g., Moffat [1978], Krause & Radler [1980]). The 
proportionality M(k) a k' V2 is indicated for the inertial sub-range of three-dimensional, 
homogeneous, isotropic, incompressible hydromagnetic turbulence; however, kinetic helicity 
injection at large k leads to an inverse cascade of magnetic helicity and M(k) a k A at low k 
[Pouquet, Frisch & Leorat, 1976; Stevenson, 1983]. The latter shows similarity between an n' 1 
spectrum and a downwardly continued observational spectrum i? n (0.55a, 1965) for n < 8. 
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Next define horizontal wave-number s [n(n+l)/c 2 ] 1/2 via the surface Laplacian. Granting 
M(k) k' 1 at low k, if M(k) « {R n c (c)} and if A: 2 k^ 2 , then {R n c (cJ} [n(n+l)]' 1/2 at low n and 

{R n c faj} = K\ [w(n+l)]" 1/2 (c/a) 2n+4 (Ala) 

= K m (n+my m (c/af D+A (Alb) 

= K s n 1 (c/a ) 2n+4 . (Ale) 

Here K\ , Km and ATs denote constants, c core radius, and a Earth’s radius. By inequality (3a) for 

n > Ad, if either Stevenson’s [1983] relation (Ale), McLeod’s [1996] rule (Alb), or equation 
(Ala) holds at low degrees n < As, then Ae < Ad, then there may indeed be an intermediate, if not 
inertial, sub-range between Ae and Ad- 

Of course, the “ifs” strain a comparison already made difficult by possible effects of mantle 
heterogeneity near the CMB; anisotropy imposed by rotation at planetary angular velocity £2, the 
CMB, and the field itself; compression; suppression of turbulence by rotation or a strong field; 
and sources in the mantle and crust. Yet as argued in Appendix B, the net effect might amount to 
a time-averaged spectrum similar to (Ala), perhaps with power distributed unevenly among the 
various orders within each R n , notably R\, due to anisotropy and lateral heterogeneity. There is 
evidence that spectra like (Ala) describe both modem observational spectra and time-averaged 
paleo-field behavior [ Voorhies & Conrad, 1996]. Moreover, agreement between such spectra 
and observation [ Voorhies et al., 2002] offers some support for the temporal magnetic power 
spectrum used by McLeod [1996] to obtain (Alb) for n > 2. 

Appendix B: Spectra from Scale Analyses 

A reviewer asks if the spectrum of a magnetic field in a turbulent fluid has physical 
plausibility analogous to the famous Kolmogorov k' scaling for the kinetic energy density 
spectrum E{k) in the inertial sub-range of a turbulent flow (see, e.g., Tennekes & Lumley [1972]) 
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and requests mathematical illustration. Here we show the answer is yes and offer an illustration 
that yields both expected spectrum (3b) and SV spectrum (5). As noted by Pouquet et al. [1976], 
however, care is needed to obtain a magnetic spectrum excited by small scale flow, instead of an 
Alfven-wave spectrum excited in an inertial sub-range maintained by large scale flow (see B.l 
below). Additional care is needed to obtain a spectrum for a self-governing dynamo, rather than 
one driven by assumed rates of kinetic energy injection (see B.2). Finally, to deduce theoretical 
spectra comparable with observational spectra, extra care is needed to account for sphericity and 
the anisotropic dynamical effects of rotation, Lorentz forces, and the CMB (see B.3 and B.4). 

Consider flow of a magnetized Newtonian fluid with scalar mass density p, kinematic shear 
viscosity v, magnetic permeability p, electric conductivity a and magnetic diffusivity 1/pa = rp 
Quasi-steady changes in macroscopic B due to fluid velocity u and magnetic diffusion are 
described by the induction equation (4). For a characteristic length scale L and flow speed U, the 
ratio of motional to diffusive terms scales as the magnetic Reynolds number 

Rm m ULJx\ = | Vx(uxB) | IqVxVxBl' 1 . (Bl) 

This is also the ratio of the characteristic eddy diffusivity to the magnetic diffusivity. 

Much as the kinetic energy transport equation is obtained from the inner product of u with 
the momentum equation, the magnetic energy transport equation is obtained via the inner product 
of B with the induction equation (see, e.g., Chandrasekhar [1981]; Gubhins & Roberts [1987]). 
The Ohmic dissipation of magnetic energy per unit volume, -J /a, scales as r\B /2pL . The 
viscous dissipation of kinetic energy per unit volume, -pvu« VXVXu for solenoidal flow, scales 
as vpU 2 /2L 2 . The ratio of magnetic to viscous dissipation scales as (v/q )"* (B 2 / pp H 2 ), where v/r) 
is the magnetic Prandtl number. In the flow of a fluid metal with v/r| « 1, if magnetic energy 
density B 2 /2p is at least as great as kinetic energy density pu 2 /2, then magnetic dissipation will be 
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very much greater than viscous dissipation. So we focus on magnetic dissipation and dynamo 
action, whereby motion of the fluid conductor across the field does work against the Lorentz 
force, converting kinetic into magnetic energy. The latter may accumulate, dissipate into heat, or 
radiate away. The rate of work done against the Lorentz force per unit volume, -u»(JXB), scales 
as UB z /2\xL. 

More generally, u and B may vary on all possible length scales l and may be represented 
mathematically via superposition of orthogonal modes, such as Fourier transforms with wave- 
vector k and wave-number I k I = k ~ 1 l . Work done by flow at one scale against the non-linear 
Lorentz force can energize the field over a range of scales. And magnetic energy at one scale is 
influenced by the flow over a range of scales. Because of this mode mixing, magnetic energy 
dissipated at any single scale may come from kinetic energy distributed over a range of scales. 

Mode mixing is governed by the usual selection rales. We shall, however, only consider 
states of a hydromagnetic system with a standard spectral deviation <[M(k) - <M(k )> ] > that 
does not vastly exceed the mean <M(k)>. In such a state, a theoretical expectation spectrum 
{M(k)} ~ <M(k)> could be of some use; moreover, there is a fair chance that a sample spectrum 
at a single time is within a factor of two or three of the mean. In or near such a statistically 
steady state, energy mixed from one mode k a to other modes k, at one rate is re-mixed to still 
other modes k s and, on average, returns to mode k a at about the same rate. This implies the 
transformed magnetic and velocity fields are so thoroughly intertwined that fields of like wave- 
number are physically related, as can be seen by repeated or multiple applications of the selection 
rales. We focus on the physical relations, with the understanding that they result from mode 
mixing rather than a single application of the selection rales. 
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Suppose u and B can be considered smooth on very small spatial scales, with scale lengths < 
X 0 . At length scale A 0 , let Rm attain unity at characteristic speed u = o 0 : 

t)(yWT| = 1 • (B2) 

At this scale and speed, eddy diffusivity ul = o 0 X 0 will equal magnetic diffusivity r\ and advective 
time-scale x a = II u will equal diffusive time-scale x d = l 2 /r\. At larger scales, or at greater speeds, 
fluid motion may curl the field faster than it diffuses away via electrical resistance to its source 
current of density J (VXB = pj, VXVXB = -V 2 B). 

Let b(l ) denote the magnetic field of length scale l and let bo denote b(k o). At / = and 
speed Do, the magnetic dissipation scales as Do(Ao) = y\bo‘l2\xko. By equation (B2), this also 
scales as the rate of work against the Lorentz force: Do ~ x\bol2\xko' = v>obo‘!2\xko. So A-o is the 
magnetic dissipation scale, \q ~ (x\bo H\xDo) , and t)o is about [2Do/ctbo ] • 

B.l Kolmogorov-Alfven Scaling 

Denote by £ the total magnetic dissipation per unit volume in a hydromagnetic flow. Suppose 
this occurs mainly over a range of length scales less than or approximately equal to ~ko. Further 
suppose the magnetic energy dissipated comes from the kinetic energy p?7 / 2 of a flow with 
speed U and large geometric length scale L » X 0 - The energy is supplied by work done against 
the field B of length scale L. If this energy conversion occurs in advective time-scale LIU , then 
e - UB 2 I2\xL - pL 3 /2L . (B3) 

More generally, the energy dissipated may come though a range of scales. Consider a sub- 
range Xq<1 <L and x>o<u<U through which magnetic energy cascades from the field of scale L 
toward the dissipation range at rate 

£ = ub 2 l2[d = p« 3 /2Z . (B4a) 
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This assumes the time-scale for kinetic energy density p u / 2 in motions u(l) of scale l to be 
transferred to the field b(l) is x a . Granting (B4a) for now and solving for u gives: 

u » (2| xleJb 2 ) * (2Ze/p) 1/3 . (B4b) 

The implied kinetic energy per unit volume at wave-number k ~ l~ l is pu 2 /2 ~ (p/2) 1/3 (e !k) m \ 
therefore, the kinetic energy density per wave-number is 

£ka(£) = pu 2 /2k * (pe 2 /2 ) 1/3 fc' 5/3 . (B5) 

This is proportional to the Kolmogorov form for an inertial sub-range. Also by relation (B4b), 
the magnetic energy per unit volume at k is b 2 /2p ~ (p/2) 1/3 (e/fc) 2/3 ; therefore, in this case the 
magnetic energy density per wave-number is 

Mka(® = b 2 l2\xk = (pe 2 /2 ) 1/3 fc' 5/3 . (B6) 

As appropriate to a field of Alfven waves excited by motions in an inertial sub-range, £ka(£) and 
Mka(T) are approximately equal. Of course, relations (B4a) though (B6) fail when dynamical 
constraints render advective time-scales x a irrelevant to the transfer of kinetic energy. 

B.2 Large Scale Flow 

By presuming the rate of work done against the field equals the kinetic energy density per 
advective time-scale x a , case B.l compels magnetic and kinetic energies to be in approximate 
balance, as in an Alfven wave field. For dynamo action, however, magnetic energy density 
b 2 (l)l 2p need only be replenished over free decay time x d . The rate of work done against the 
Lorentz force still scales as ub 2 /2pl, but the kinetic energy transfer only needs to occur over time 
l /r). In this case, relation (B4a) would be replaced with 

8 = ub 2 /2pl * pw 2 r|/2Z 2 . (B7a) 

Granting this for now and solving for u yields 

u » (2/ 2 E/pri) 1/2 = (2 pl^b 2 ) = (lb 2 ipr\\x) . (B7b) 
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This relation implies 

b 2 /2ii « (pr|e/2) 1/2 . • (B8) 

By relation (B8), the magnetic energy density at k ~ T 1 is also about (pr|£/2) 1/2 ; therefore, in this 
second case, the magnetic energy density per wave-number is 

Mlsf(£) = b 2 ll\ik * (py\eJ2) m k'\ (B9) 

Also from relation (B7b), the kinetic energy density at k is p« 2 /2 = (dr\)k' 2 \ therefore, in this 
case the kinetic energy per unit mass per wave-number is 

£lsfW = P« 2 /2fc * (s/r|)^' 3 . (BIO) 

This k~ 3 form describes a predominantly large scale flow. It illustrates how a strong field may 
suppress small scale flow. Indeed, E(k)/M(k) from relations (BIO) and (B9) is (2£/pr| 3 ) 1/2 £' 2 , 
which falls off as k~ 2 and increases with o 3/2 and £ 1/2 . The spectrum of magnetic change, obtained 
using relations (B9) and (BIO), the motional term in equation (4), and implicit treatment of mode 
mixing, is proportional to that of the field itself: Fi^?(k) = (d t b) 2 /k ~ ( kub) 2 /k ~ (2£p) 3/2 (a/p) 1/ V 1 . 
Evidently, this is not the case, near the top of Earth’s core. 

Case B.2 might seem to describe a more efficient dynamo than case B.l because, presuming 
Td » x a , the kinetic energy required to maintain a large scale field would be extracted over a 
longer time. In fact, case B.2 still presumes: (i) Ohmic dissipation is confined to small scales < 
h) on which Td < x a ; (ii) an energy cascade from large to small scales; and (iii) a time-scale for 
kinetic energy transfer devoid of dynamical justification. 

B.3 Small Scale Flow 

Now suppose most of the kinetic energy is injected at fairly small scales near l*, albeit with 
l* > Xq. This motion may set up a cascade of magnetic energy to smaller scales, but may also 
drive a mode-mixed reverse cascade to larger scales. Let £ indicate the overall magnetic 
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dissipation in the flow, which realistically occurs over the full range of length scales. At larger 
scales l > l*, and near a statistically steady state, the rate of work done against the Lorentz force 
approximately balances Ohmic dissipation, so 
|u*(JXB)| = a'u(l > > l*)] 2 /\xl 

» [b(l > l*)f/\i 2 al 2 - I J 2 /a I . (B 1 la) 

Here a is a constant representing typical geometric factors; it is large if the field tends to be 
nearly Lorentz force-free, or if u tends to be parallel to either J or B. We solve (Blla) for 
u(l > l*) ~ ar\/l , (Bllb) 

and see a as a pseudo-scale invariant magnetic Reynolds number. By (Bllb), we expect 

{E(k <k*)} = (p[w(£ < k*)f/2k} * a 2 pr| 2 &/2 . (B12) 

This is consistent with a predominantly small scale flow. 

We cannot assume the time-scale on which work is done at the expense of kinetic energy at 
length scale l because it depends on dynamical constraints near the top of a self-excited, self- 
regulating core geodynamo. Following Voorhies [1991, equation (B5)] and Benton [1992, 
equation (32)], the dynamical constraints are here summarized by the magneto-geostrophic radial 
vorticity balance near the top of the core, 

pVs*(2£2 r u s ) = -Vs*(B r J s ), (B13) 

where Vs* denotes the surface divergence operator, Q r the radial component of planetary angular 
velocity (£2 o cos0), B r the radial component of B, and u s and J s the horizontal components of u 
and J, respectively. Anisotropies imposed by rotation and by the field are explicit in equation 
(B13). The published derivations explicitly account for anisotropy imposed by the CMB by 
omitting small terms with u r and J r , terms which would be zero at the spherical boundary with a 
rigid, insulating mantle. The condition J r = 0 implies any poloidal current, hence toroidal 
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magnetic field, is zero at the boundary; therefore, a magnetic spectrum deduced using (B13) 
1 should only be applied to the poloidal field at the top of the core and, upon upward continuation, 

to the potential field of core origin. 

If the simple scaling of (B13) were multiplied by l, it would suggest Lorentz and Coriolis 
forces are of similar magnitude. Yet (B13) concerns vertical vorticity, not forces, and can hold 
when the Lorentz force is weak. Moreover, such scaling overlooks some geometric effects. For 
example, axisymmetric zonal flows contribute nothing to the left side of (B13), yet may be 
important in the force balance; similarly, a class of currents tangent to contours of B T may help 
generate the main field, yet contribute nothing to the right side of (B13) [Voorhies, 1991]. To 

i 

avoid any misimpression, equation (B13) is here scaled as 

2pp £2u(/ >/*)// = [b(l > l*)] z /\xl 2 , (B14) 

where (3 is a constant representing typical geometric factors. On the hypothesis of a very weak 

i 

Lorentz force near the top of the core, (3 « 1 . 

From relation (B 14), we find 

[b(l > Z*)] 2 /p * 2(3p Qlu(l > l *) . (B15) 

We use relation (B 1 lb) to eliminate u from (B 1 5) and obtain 

l 

[b(l>l*)] 2 /\x * 2appQri . (B16) 

i 

From this relation, the expected magnetic energy density per wave-number is 

{M(k<k*)} = {{b(k<k*)] 2 /2pk} * aPpGri** 1 . (B17) 

If relation (Blla) between Ohmic dissipation and the rate of work done against Lorentz 
| forces can be extended to a rate of change of kinetic energy, so that r\b 2 /pl 2 ~ ub 2 /\xl ~ pu 2 /x, we 

can solve for the transfer time-scale x = (pu 2 \i/b 2 )l 2 /r\. By relations (Bllb) and (B1 6), this 
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dynamical time-scale turns out to be t = a/2Qp, another pseudo-scale invariant for l > l *. Of 
course, for a » 1 and 3 « 1, t » 

With the right hand side of induction equation (4) scaled according to the motional term, and 
using relations (B12) and (B16), the expected SV spectrum is 

{F(k >k*)} = {0 t bflk) = [ku(k < k*)b(k < k*)] 2 /k = 2a 3 3ri 3 ppQ k?. (B18) 

The ratios {2 \xM}/{F} from relations (B18) and (B17) define scale-variant squared time-scales 
(ar|& 2 )" 2 which constrain a independent of 3 (see section 6). 

If l* ~ Xo, an energy cascade from k* to larger k would go directly into the dissipation range. 
Instead suppose l* > Xo and consider the intermediate sub-range Xq < l < l*. Let e s denote the 
magnetic dissipation on scales < Xo> so e s < 8. This dissipated energy is re-supplied by work done 
against the Lorentz force which, in turn, comes from kinetic energy on scales near l*. This 
energy input is now at relatively larger scale l* > l, so one might re-consider cases (B.l) or (B.2). 
Motions in the sub-range should contain considerable kinetic energy, perhaps suggesting use of 
(B4a). By relation (B17), however, these motions are embedded in a relatively larger scale field, 
which may favor (B7a). Yet a small scale eddy cannot be everywhere parallel to a large scale 
field, so neither the rate of work against the Lorentz force nor the magnetic dissipation will be 
negligible in the intermediate range. Therefore, neither case (B.l) nor (B.2) need apply. Further 
analysis of an intermediate, if not inertial, sub-range is omitted for brevity. 

B.4 Simple Discretization and Accounting for Sphericity 

We need to relate continuous Cartesian magnetic spectrum (B17) with discrete spherical 
harmonic spectra (Ala), (3c) or (11c), and SV spectrum (B18) with (5), albeit only to the order of 
magnitude accuracy appropriate to a simple scale analysis. To do so, first integrate continuous 
spectrum (B 17) over the small domain from k - AkJ2 to k + Ak/2, 
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(B19a) 


k + Ak/2 k + Ak/2 

f {M(k> k 0 )}dk - a(3pQri \ k A dk 

k - Ak/2 k - Ak/2 

= app&ri ln[( 1 + Ak/2k)/(l - Ak/2k)] (B19b) 

= app£2ri(A kJk) . (B19c) 

where the logarithm has been expanded and approximated in terms of Akl2k « 1. 

Next, guided by the surface Laplacian operators in the transform domains, set k = n(n+l)/c . 
This implies A k = (n +l/2)[n(n + l)]' I/2 c' l An, and we set An = 1 to express the magnetic energy 
density in harmonics of integer degree n. With these identities we obtain 

k + Ak/2 

j {M(k < ko)}dk * a(3p£2r| (n +l/2)[n(n + l)]' 1 . (B20) 

k - Ak/2 

Note k has been treated as a horizontal wave-number k^. Distinctions between radial wave- 
vector component k r and k^ are omitted for three reasons. First, we are concerned with the 
spectrum of the field on a thin shell by the top of the core. Second, apart from the tiny jump in 
horizontal components across a viscous sub-layer, this field should match a core-source potential 
field (solenoidal B = -VV, so -V 2 V = 0, k T 2 + fc h 2 = 0, and k r = ±ik h ). Third, the vertical vs. 
horizontal anisotropy implied by CMB conditions was already used to derive equation (B 13), 
hence is implicit in relation (B 17). This omission of distinctions between kr and k^ is consistent 
with the constant aspect ratio hypothesis (see section 6). 

Relation (B20) gives the magnetic energy density per harmonic degree. The expected 
squared magnetic field per harmonic degree at the top of the core is 2\x times this quantity, so 
{R n (c')} - (2a3)(pf2pr|) (n +l/2)[n(n + l)]' 1 (B21a) 

at low degrees. On upward continuation to radius r > c, this yields 


{/? n c (r)} = K(« + ‘/2)[»(n+l)]- 1 (c7r) 2n+4 , 


(B21b) 



which is relation (3b). Though obtained in a very different way, it is indistinguishable from 
equation (lib) and differs but slightly from equations (Ala), (Alb), and (Ale). We have, 
however, identified K = 2a|3pp,r|£2. Granting pi = |0o, p = 10 4 kg/m' 3 , Q = 7.29x1 CT 5 s' 1 , and r| = 
1.6 m 2 /s, the estimate K = 4.5xl0 10 nT 2 [Voorhies et al., 2002] implies a(3 is roughly 1.5xl0' 2 . 
With a~R m ~92 from Appendix G, we find (3 = 1.7xl0’ 4 « 1. This is consistent with, and 
arguably requires, a dynamically weak field near the top of the core. These values further yield a 
dynamical time scale ct/2Q(3 of order 120 years. 

Similar operations on relation (B 18), albeit with no multiplication by 2\x, yields 
{F n (c)} * 4a 3 (3ri 3 ppQ c 4 (n + U2)[n(n + 1)] (B22a) 

for low degrees, or, on upward continuation to r > c, 

{ F-n(r ) } = 4a 3 (3q 3 ppQ c*n(n + 1/2) [n + 1] (c'/r) 2n+4 . (B22b) 

This is indistinguishable from equation (5). The ways in which these relations are obtained could 
hardly be more different, but both suppose narrow scale flow and so give the same result. 

Appendix C: Derivation of an Expectation SV Spectrum 
The magnetic change induced by each narrow-scale eddy can be approximated in the far field 
by an equivalent source of change atop a mainstream of radius c. Though a single magnetic flux 
vector on c acts as the point source of an offset dipole field, lateral advection replaces it with an 
adjacent vector of slightly different orientation and magnitude in time At; therefore, the net 
change in the exterior field is equivalent to that from a differential change in magnetic dipole 
moment Ad. Here we derive the secular variation (SV) spectrum from one such equivalent SV 
source, then the SV spectrum (5) expected from I random, uncorrelated equivalent SV sources. 

Denote by AVj(r) the change in scalar magnetic potential at observation point r due to the 
equivalent magnetic dipole moment change Adj at position rj corresponding to the i th eddy. Note 
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I r I > I r;| = c and index i = 1, 2, 3, . . I. Further denote by D l the normalized rate of change 
Adj/jioAt. The superposition of all I such changes gives the total change in potential at r, 

I ID' 1 

AV(r) = Z AVi(r) = -Z — *V At (Cl a) 

z=l i=l 4n ! r - r, I 

I 

= Z (DV4tt) • V* | r - r, I ' : At (Clb) 

z=l 

where V* denotes the gradient operator in terms of, and acting on, r; coordinates (radius q, co- 
latitude 9i, east longitude <j)j). With the Schmidt normalized associated Legendre function of 
degree n and order m denoted P n m (cos0), the spherical harmonic expansions of the AVj are used 
to rewrite (Clb) as 

/ oo n 

AL(r) = Z aAt Z (a/r) n+1 Z [A n m i cosm<\> + 5 n m , sinm<t>] P n m (cos0) . (Clc) 

i=l zz=l m= 0 


To determine coefficients (A n m j, B n m i) in expansion (Clc) from moment changes D\ recall 
°° n (n) n 

I r - ril' 1 = Z Z [C n m (0,<l))C n m (0 i ,<j) i ) + S n m (0^)S n m (0i^i)] , (C2) 

n=l m - 0 r D+1 

where C n m = cosm(j)P n m (cos0) and 5 n m = sinm4>f , n m (cos0) (see, e.g., Jackson [1975, eqn. (3.70) or 
Langel [1987, eqn. (195)]). With the components of D‘ denoted (D r ‘, Dq, £><p’), we use (C2) to 
evaluate the gradient in (Clb), equate the result with (Clc), and obtain 

a 

A n m i = (47ta 3 )‘ 1 (r i /a) n - 1 [A i nC n m (0 i ,(|) ,) + Dj — C n m (0 \,$) - A;(m/srn0 1 )5 n m (0 1 , (}),)] (C3a) 

39i 

3 

B n m i = (47ta 3 )' 1 (r, 7 a)"' 1 [D r i n5 n m (0 i ,(j) i ) + Dj —S u m {Q j.fo) + Dj (m/smOOC^O;,^)] . (C3b) 

a0j 
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In the limit as At approaches zero, the equivalent core-source SV coefficients are 


I 

d t gn m = s A n ffl i (C4a) 

i=l 

I 

d t h n m = 2 BrTi . (C4b) 

i=l 

Substitution of these coefficients into equation (2) gives the equivalent core-source SV spectrum. 
Clearly, statistical hypotheses about moment changes Ad;, notably those implied by physical 
hypotheses about eddy induced SV, can be used to derive an expectation core-source SV spectra 

n I I 

{F n c (»n;} = (n + l)(a/r) 2n+4 { 2 [( 2 An™) 2 + ( 2 B^f) (C5a) 

m= 0 i=l i=l 

nil 

{F n c (r>n)} = (n + l)(a/r) 2n+4 { 2 [( 2 2 (An m iA n m k ) + (B n m , B n \)} } (C5b) 

m- 0 i=l k=l 

Suitably correlated, perhaps non-conserved, equivalent SV sources may also describe effects 
of molecular magnetic diffusion; unfortunately, the linearity of the diffusion operator might lead 
one to an equivalent source representation of the core field itself, which is judged ill-advised. 

SV Spectrum from a Single Change in Dipole Moment 
The SV spectrum from a single D‘ is 
n 

Fa(r > n) = (n + 1 )(a/r) 2n+4 2 (A n m i) 2 + (B^) 2 . (C6) 

m= 0 

By (C3a) and (C3b), the sum over order m in (C6) is 

n n d 

2 (A n m 0 2 + (B n m i) 2 = G n 2 [(D r 'nP a m ) 2 + (Dj — P n m ) 2 
m = 0 m - 0 d0, 

d 

+ (Djm P n m /sinQi) 2 + D^Dq n P„ m (~ P a m ) 2 ] (C7) 

d0i 
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where G n (r0 = (4na 3 y 2 (r, fa) 2n ' z and both P n m and its derivative are evaluated at cosGi. 


To cast (C7) into a more illuminating form, we use spherical harmonic identities to reduce 
the four terms on the right, each a sum over m. For arbitrary (0,<)>), the sum rule for Schmidt 
normalized harmonics (see, e.g., Jackson [1975, eqn (3.62) with y = 0], Langel [1987]), 


2 [C n m (0,<t>)] 2 + [S n m (0,<j))] 2 = Z [P n m (cosQ)] 2 = 1 , 

m - 0 m-0 


(C8a) 


is differentiated repeatedly with respect to 0 to obtain 


n d 

Z 2 P n m —P n m = 0 ; 
m - 0 d0 


(C8b) 


n d 2 d 

Z 2P n m P„ m + 2[ — P n m ] 2 = 0 . 

m=0 d 2 0 d0 


(C8c) 


The first term on the right of (C7), the radial term, sums to (P/n) 2 by (C8a); the fourth term, the 
cross term, sums to zero by (C8c). 

To reduce the second, co-latitudinal term on the right of (C7), recall that r 2 times the surface 


Laplacian of S n m gives 

V s 2 S n m = -n(n+l)S a m = (sin0)‘'9 e [sin0 3 9 5 n m ] + (sin0)- 2 3 (p 2 5 n m ]. 
By adding this to the corresponding relation for C n m , we find 


(C9a) 


d 2 cosG d 

P n m = [-n(n+l) + (m/sin0) 2 ]P n m P n m ; 

d 2 0 sin0 d0 


(C9b) 


therefore, 


n d n cos0 d 

S P n m —Pn m = s [(m/sin0) 2 - n(n + l)][P n ra ] 2 P n m --P n m • 

m = 0 d0 m = 0 sin9 d0 


(C9c) 


The last term on the right sums to zero by (C8b); substitution of the remainder into (C8c) yields 
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n d 

2 [— P n m ] 2 

m - 0 d0 


(CIO) 


n 

= I [n(n + 1) - (m/sin0) 2 ](P n m ) 2 . 
m = 0 

To reduce the third term on the right of (Cl) and complete reduction of the second, we need 
the sum over m of -[mP n m (cos0i)/sin0i] 2 In terms of the angle y between the r and r; , x = cosy, 
and Legendre polynomials P R (x), it has been shown that this sum reduces to 


n dP n (x) 

2 -[mP n m (cos0i)/sin0i] 2 = | x = ) . (Clla) 

m - 0 dx 


Evaluation of Legendre’s equation at a: =1 and the normalization P n (l) = 1 then imply 

n 

2 -[mP n m (cos0i)/sin0i] 2 = - n(n + l)/2 . (Cl lb) 

m = 0 

The proof of (Cl lb) by Whaler & Gubbins [1981, Appendix B], which they attribute to P. H. 
Roberts, is more elegant than that by Voorhies [1998]. 

By (Cl lb), (CIO) is equal to n(n+l)/2; therefore, the second term on the right of (Cl) sums to 
n(n + l)(De') 2 /2. Similarly, the third term on the right of (Cl) sums to n(n + l)(D^,') 2 /2. 
Substitution of (Cl), as simplified by identities (C8b,c), (CIO) and (Cl lb), into (C6) yields 

F n c (r > n) = (47ir 3 )- 2 (r,/r) 2n - 2 (n 2 [n + l](D r ‘) 2 + n[n + l] 2 [(D e i ) 2 + (D* 1 ) 2 ]) . (C12) 

This is the SV spectrum from a single changing moment Adj/At = poD 1 at n < r. 

Expected SV Spectrum from Uncorrelated, Randomly Oriented Changes in Moment 

At an instant in geologic time, the I equivalent SV sources D 1 representing eddy induced SV 
are treated as a sample population with mean square moment change {D 2 }. The coefficients in 
(C5a,b) are given in terms of the D 1 by (C3a,b). The source positions (0i,<t>O on c are taken to be 
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random samples from a laterally uniform spatial distribution. Such positions are expected to be 
uncorrelated in that {ri»r k } = 8 ik {(r0 2 }, where the Krpnecker 5 ik is 1 if i-k and is 0 if i * k. 

For randomly oriented moment changes corresponding to a kinematically unbiased ensemble 
of eddies, any particular orientation is as likely as its opposite, 

{D r '} = {ZV} = {Dty} = 0, (C13a) 

there is no reason to expect cross-correlated components for an individual moment change, 
{AW} = {DjDj} = {DqD$} = 0, (Cl 3b) 

there is no reason to expect cross-correlated moments 

{ZVA k } = {A ! A k } = { DjD^} = 0 for i*k (C13c) 

{D r 'De k } = {D t 'D^} = {DqD^} = 0 for i * k, (Cl 3d) 

but the auto-correlations remain perfect 

{D;d x 1 } = {DjDj} = (W) = {D 2 }/3. (C13e) 

Equations (C13a-e), summarized by {D/} = 0 and {Dj‘Di k } = { D 2 /3 } 5i k 5ji, provide a 
mathematical statement of “random equivalent SV sources” associated with a kinematically 
unbiased ensemble of eddies. Polarizing dynamical effects, whereby eddies might tend to align 
with the planetary rotation axis and/or the field, are outside the focus of this appendix. 

Though D 1 depends on position (0i,4>i), position is independent of D 1 ; therefore, when 
evaluating {A n m i} and {B n m i} via (C3a,b), the expectation operator passes through the harmonic 
functions of (6i,(j)i). These expectation values are zero by (C13a), as are expected or mean SV 
coefficients by {(C4a,b)}, hence the expected components of SV itself. Though the expected SV 
vector from random equivalent sources is zero, the expected spectrum is not. 

Equations (C13a-e) and (C3a,b) imply contributing coefficients are not cross-correlated 
{An” AA} = { A” i Bn m k} = o for i*k. (Cl 4) 
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Again, the auto-correlations remain. Substitution of (C14) into (C5b) yields 

n I 

{F n c (r > r\)} = {n + l)(a/r) 2n+4 { Z [ £ (A n m ,) 2 + (B^i) 2 ] } (C15) 

m-0 i=l 

The order of sums in (Cl 5) is reversible, so the expected spectrum from I uncorrelated SV 
sources is the sum of expected spectra from each source. With r\ = c, (C15) and (C12) imply 
{F n c fr > c)} = (47tr 3 )' 2 / (c7 r) 2n ’ 2 n(n + l)[n{D r 2 } + (n+ 1)({D 9 2 }+ {De* 2 })] (C16) 

By (C13c), this simplifies to 

{F n c (r > c)} = (4nr 3 )' 2 (I{D 2 }/3) (c7r) 2n ‘ 2 [n(n + 1)(2 n + 1)] (C17) 

With C = (47tc' 3 )’ 2 (2/{£> 2 }/3), (Cl 7) becomes 

{F n c (r > c)} = Cn(n + Vi){n + 1 )(cVr) 2n+4 , (C18) 

which is just equation (5). 

Appendix D: SV Spectral Variance 

For steady c, equation (5) also gives the time-averaged spectrum expected from a temporal 
sequence of kinematically unbiased ensembles of equivalent SV sources, denoted <{F n c (f)}>, 
albeit with amplitude C replaced by <C(t)>. The associated sequence variance <[{F n c (r)} - 
<{F n c (X)}>] 2 > is proportional to <[C(t) - <C(t)>] 2 >, but has the same functional dependence on n 
as the square of spectrum (5). The root mean square sequence deviation is therefore proportional 
to the expected SV spectrum itself. This is a strong restriction on the kinds of long term 
magnetic variations that could be adequately described by the narrow scale transport model. For 
example, if eddies and eddy transport change in response to a change in a driving function (e.g., 
heat flux across a CMB, heat and buoyant component fluxes across an inner core boundary, etc.), 
then this sequence of kinematically unbiased ensembles could describe a corresponding change 
in the amplitude, but not the shape, of the expected low degree SV spectrum. 
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Different sets of Mi can represent broad-scale SV induced by different numbers of eddies 
transporting different fields at different speeds; however, a single ensemble from this temporal 
sequence of ensembles will give but one ensemble-mean low degree SV spectrum. Each such 
ensemble represents lateral magnetic transport by an ensemble of eddies with a corresponding 
mean number of eddies, mean speed and, arguably, a mean pattern and intensity for the field 
being transported. If the mean number of eddies, mean speed, and mean field properties can and 
do change by their own magnitude from one ensemble to the next in the sequence, say from one 
epoch to another, then the sequence variance will not only be proportional to, but will arguably 
have amplitude similar to, or perhaps somewhat greater than, <{F n c }> 2 itself. 

Physical deviations F n (t) - {F n c } may also arise from dynamic processes that cause serial, as 
well as lateral, correlation between SV induced by eddies which grow, vacillate and decay. Such 
deviations are arguably of a magnitude similar to {F n c } itself and may persist for a long time, 
perhaps as long as dynamic JxB conditions associated with a polarity chron, a particular pattern 
of heterogeneous electrical or thermal boundary conditions, etc.. If so, and if spectrum (5) still 
described a geologically longer time average of F n °(t), then natural fluctuations F u c (t) - {F n c } 
would again have variance similar to, if not somewhat greater than, {F n c } 2 . 

Appendix E: Effects of Uncertain Process Variance 
The statistical significance of residuals depends on both data covariance and process 
covariance. Here, data covariance represents uncertainties in an observational SV spectrum 
F a (a,t), hence the sparse distribution of uncertain geomagnetic measurements used to determine 
SV coefficients. Process covariance represents fluctuations of a theoretical core-source SV 
spectrum F n c (a,t) about its expectation value, plus omitted non-core contributions to the 
observations. For F n (a,t) determined by careful analyses of satellite and surface observatory 
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measurements, and anticipating process variance of magnitude similar to {F n c (a )} 2 itself, process 
variance should exceed data variance at low degrees. 

The narrow scale flow hypothesis specifies neither a probability distribution function for SV 
spectrum fluctuations nor a process variance, so one can but estimate the latter from residuals. I 
doubt this distribution can be determined from 10 or 12 sample residuals. Fluctuations F n c (t) - 
{F a c } cannot be normally distributed because F n c (t) is positive, but many other distributions are 
eligible (truncated Gaussian, log-normal, chi-squared, etc.). Whatever the trial distribution, when 
process variance is estimated from but a single set of residuals, here by scaling weights so the 
sum of square weighted residuals per degree of freedom becomes unity, then the (in)significance 
of the residuals does not test the hypothesis. The test is in the comparison of estimated 
parameters with values determined by independent analyses of independent data. We therefore 
emphasize comparison of magneto-spectral estimates of source radius c with independent c s . 

Appendix F: A Softer SV Spectrum for 2000? 

On orbit magnetic calibration and magneto-optical alignment checks enabled the fitting of 
high quality main field and SV coefficients to high precision scalar and oriented vector data 
acquired by 0rsted. Here we consider F n (2000) from model OSVM [Olsen, 2002], which fitted 
24,585 two-component vector and 68,448 scalar and field aligned 0rsted data from a 2.5 year 
interval and data from about 100 observatories for 1998-2000. Compared with <F n > for 1960- 
80, there are decreases in F 2 and perhaps F 6 , increases in Fj, F 4 , and F 5 , but rather small changes 
for Fi, F%, and F 9 . The decrease in F 2 suggests less rapid quadrupole rebound. 

Estimates of c using OSVM and spectrum (5) are typically well below c s . Estimates using 
proxy spectrum (7) and rc„ „ n = 1, shown in Tables Fla, agree fairly well with c s and have less 
misfit than shown in Tables 3a and 5a; moreover, for = 3, the errors shown in Table Fib are 
rather less than found in Tables 3b and 5b. We therefore suggest spectrum (7) might be more 
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suitable than (5) at some epochs, perhaps due to changes in lateral length scales of core surface 
flow, hence in the lateral correlation length scales of SV at the top of the core. Whether or not 
this is substantiated by other epoch 2000 models with non-steady SV, different filters of 
observatory data, normal weights, or different treatments of external fields remains to be seen. 

A broadening in scale of core surface flow between the interval 1960-1980 and the 0rsted 
epoch 2000 may explain a softer SV spectrum, though other explanations are possible. If 
spectrum (5) were be shown to be wholly inadequate outside the 1960-1980 interval, yet (lib) 
remained adequate, then arguments leading to (lib) might have to be revised - possibly by 
including effects of Lorentz forces on eddy diffusivities. In this way it might be possible to 
discern deviations from tangential geostrophy by spectral methods. 

Appendix G: A Bound on Core Conductivity 

Until a few years ago, observational SV spectra were not thought to be well enough 
determined to high enough degrees to test spectrum (5) against seismology; it was only checked 
via (10b) in the context of (lib) [Voorhies & Conrad, 1996]. From the check of (10b), we infer 
there is but a single important lateral empirical diffusivity, independent of degree n. Moreover, 
from the value of oto, this lateral empirical diffusivity is £ = c 2 /oq = 147 m 2 /s with an uncertainty 
factor or divisor of two. If this is equal to a single ensemble mean lateral eddy diffusivity UL , 
then there is but one important eddy magnetic Reynolds number R m = \xaUL. If one assumes a = 
5xl0 5 S/m and vacuum permeability, then this magnetic Reynolds number is R m = 92 with an 
uncertainty factor of two. Another interpretation, which does not assume a and allows forced 
dipole decay, uses oto to bound core magnetic diffusivity r\ < 2(c/7t) 2 /oco = 30 m 2 /s, and, for 
vacuum permeability, core electric conductivity a > 2.7 x 10 4 S/m, both with an uncertainty 
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factor of two. This is consistent with a fluid metallic core and, in turn, provides a check on the 
larger theory. 
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Table 1. Core radius estimated by fitting log-theoretical SV spectrum (5) to 
log-observational spectrum GSFC 9/80 for various epochs and degree ranges. 


Epoch 

n 

c (km) 

n 

1980 

1-13 

3441 ± 96 

3-11 

1975 

1-13 

3472 ±95 

3-11 

1970 

1-13 

3475 ± 99 

3-11 

1965 

1-13 

3490 ± 97 

3-11 

1960 

1-13 

3488 ± 87 

3-11 

Avg. 

1-13 

3470 ±91 

3-11 


c (km) 

n 

c (km) 


3429 ±135 

3-10 

3280 

± 127 

3539 ±106 

3-10 

3396 

± 

79 

3540 ±115 

3-10 

3387 

± 

89 

3547 ±115 

3-10 

3397 

± 

90 

3484 ± 99 

3-10 

3334 

± 

55 

3511 + 101 

3-10 

3360 

± 

61 


Table 2a: Core radius and scaled uncertainty 
from equation (5) and the time-averaged 
observational SV spectrum GSFC 9/80. 


Degrees 

q 2 

c 

Error 

Fitted 


km 

km 

1-6 

0.2316 

3026 ±174 

-454 

1-7 

0.1868 

3045 ± 124 

-435 

1-8 

0.1982 

3130 ±108 

-350 

1-9 

0.1938 

3184 ± 90 

-296 

1-10 

0.1710 

3196 ± 73 

-284 

1-11 

0.3591 

3318 ± 95 

- 162 

1-12 

0.4138 

3391 ± 91 

- 89 

1-13 

0.5041 

3470 ± 91 

- 10 


Table 2b: As Table 2a, but n m w = 3. 


Degrees 

q 2 

c' 

Error 

Fitted 


km 

km 

3-6 

0.0921 

3452 ±234 

- 28 

3-7 

0.0833 

3315 ± 151 

- 165 

3-8 

0.0709 

3378 ±108 

- 102 

3-9 

0.0579 

3398 ± 77 

- 82 

3-10 

0.0551 

3360 ± 61 

- 120 

3-11 

0.2005 

3511 ±101 

+ 31 

3-12 

0.2168 

3579 ± 92 

+ 99 

3-13 

0.2612 

3658 ± 89 

+178 
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Table 3 a: Core radius and scaled uncertainty 
from equation (7) and the time-averaged 
observational SV spectrum GSFC 9/80. 


Degrees 

q 2 

Core Radius 

Error 

Fitted 


km 

km 

1-6 

0.2040 

3488 ±188 

+ 8 

1-7 

0.1652 

3462 ±133 

- 18 

1-8 

0.1526 

35 19 ±106 

+ 39 

1-9 

0.1359 

3547 ± 84 

+ 67 

1-10 

0.1210 

3532 ± 68 

+ 52 

1-11 

0.2466 

3642 ± 86 

+ 162 

1-12 

0.2701 

3701 ± 80 

+ 221 

1-13 

0.3223 

3767 ± 79 

+ 287 

Table 3b: As Table 3a, but nmin 

= 3. 

Degrees 

q 2 

Core Radius 

Error 

Fitted 


km 

km 

3-6 

0.0791 

3826 ±240 

+ 346 

3-7 

0.0846 

3644 ± 168 

+ 164 

3-8 

0.0668 

3688 ±114 

+ 208 

3-9 

0.0534 

3687 ± 81 

+ 207 

3-10 

0.0591 

3628 ± 68 

+ 148 

3-11 

0.1743 

3774 ± 102 

+ 294 

3-12 

0.1788 

3832 ± 89 

+ 352 

3-13 

0.2070 

3902 ± 85 

+ 422 


Table 4a: Core radius and scaled uncertainty 
from equation (5) and the time-averaged 
constrained SV spectrum CM3. 

Degrees q 2 Core Radius Error 


Fitted 


km 

km 

1-6 

0.2305 

2994 ± 172 

-486 

1-7 

0.1950 

3045 ± 127 

-435 

1-8 

0.2017 

3127 ±108 

-353 

1-9 

0.2064 

3192 ± 94 

-288 

1-10 

0.1848 

3211 ± 78 

-269 

1-11 

0.1778 

3242 ± 65 

-238 

1-12 

0.1836 

3279 ± 59 

-201 

1-13 

0.3096 

3200 ± 66 

-280 
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Table 4b: As Table 4a, but ftmin — 3. 


Degrees 

Fitted 

q 2 

Core Radius 
km 

Error 

km 

3-6 

0.0743 

3408 ± 208 

- 72 

3-7 

0.0546 

3342 ± 124 

- 138 

3-8 

0.0458 

3390 ± 87 

- 90 

3-9 

0.0401 

3423+ 65 

- 57 

3-10 

0.0373 

3394 ± 51 

- 86 

3-11 

0.0321 

3400 ± 39 

- 80 

3-12 

0.0322 

3420 ± 34 

- 60 

3-13 

0.2786 

3282+ 83 

- 198 


Table 5a: Core radius and scaled uncertainty 
from equation (7) and the time-averaged 
constrained SV spectrum CM3. 


Degrees 

Fitted 

q 2 

Core Radius 
km 

Error 

km 

1-6 

0.2076 

3452 ± 188 

- 28 

1-7 

0.1664 

3462 ±133 

- 18 

1-8 

0.1517 

3516 ±106 

+ 36 

1-9 

0.1400 

3555 ± 86 

+ 75 

1-10 

0.1228 

3549 ± 69 

+ 69 

1-11 

0.1103 

3559 ± 56 

+ 79 

1-12 

0.1045 

3578 ± 48 

+ 98 

1-13 

0.3050 

3474 ± 71 

- 6 

Table 5b: As Table 5a, but n^n 

= 3. 

Degrees 

Fitted 

q 2 

Core Radius 
km 

Error 

km 

3-6 

0.0615 

3777 ± 209 

+ 297 

3-7 

0.0514 

3673 ± 132 

+ 193 

3-8 

0.0398 

3701 ± 88 

+ 221 

3-9 

0.0323 

3714 ± 63 

+ 234 

3-10 

0.0371 

3665 ± 54 

+ 185 

3-11 

0.0324 

3654 ± 42 

+ 174 

3-12 

0.0289 

3662 ± 34 

+ 182 

3-13 

0.3198 

3502 ± 94 

+ 22 
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Table 6: Core radius and scaled uncertainty 
from equation (lib) and observational main 
field spectrum OIFM 


Degrees 

s 2 

Core Radius 

Error 

Fitted 


km 


km 

1-4 

0.5588 

3441 ±575 

- 39 

1-5 

0.3730 

3461 ±334 

- 19 

1-6 

0.2803 

3444 + 218 

- 36 

1-7 

0.2326 

3496 ±159 

+ 16 

1-8 

0.2021 

3454 ±120 

- 26 

1-9 

0.2054 

3524 ±103 

+ 44 

1-10 

0.1797 

3522 ± 

82 

+ 42 

1-11 

0.1602 

3516 ± 

67 

+ 36 

1-12 

0.1539 

3542 ± 

58 

+ 61 

1-13 

0.1826 

3589 ± 

57 

+109 

1-14 

0.1782 

3610 ± 

51 

+130 

1-15 

0.3056 

3681 ± 

61 

+201 

1-16 

0.4436 

3751 ± 

68 

+271 

1-17 

0.8074 

3853 ± 

86 

+373 

1-18 

1.2241 

3958 ± 

99 

+478 


Table 7: Core radius and scaled uncertainty 
from equation (1 lb) and observational main 
field spectrum CM3. 


Degrees 

s 2 

Core Radius 

Error 

Fitted 


km 


km 

1-4 

0.7210 

3467 ±658 

- 13 

1-5 

0.4817 

3498 ±384 

+ 18 

1-6 

0.3613 

3503 ±252 

+ 23 

1-7 

0.2900 

3520 ± 179 

+ 40 

1-8 

0.2773 

3432 ±140 

- 48 

1-9 

0.2810 

3513 ±120 

+ 33 

1-10 

0.2466 

3504 ± 

96 

+ 24 

1-11 

0.2194 

3508 ± 

78 

+ 28 

1-12 

0.1975 

3507 ± 

65 

+ 27 

1-13 

0.1996 

3539 ± 

59 

+ 59 

1-14 

0.1956 

3562 ± 

52 

+ 82 

1-15 

0.2546 

3612 ± 

54 

+132 

1-16 

0.4281 

3687 ± 

65 

+207 

1-17 

0.6684 

3770 ± 

76 

+290 

1-18 

0.9795 

3859 + 

87 

+379 
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Table 8: As Table 7, but with n^n = 3. 


Degrees 

Fitted 

s 2 

Core Radius 
km 

Error 

km 

3-6 

0.0175 

3452 ±102 

- 28 

3-7 

0.0126 

3482 ± 62 

+ 2 

3-8 

0.0521 

3336 ± 91 

-144 

3-9 

0.1144 

3484 + 111 

+ 4 

3-10 

0.0958 

3476 ± 83 

- 4 

3-11 

0.0828 

3486 ± 65 

+ 6 

3-12 

0.0725 

3487 ± 52 

+ 7 


Table Fla: Core radius and scaled uncertainty 
estimate from equation (7) and observational 
SV spectrum OSVM. 

Degrees q 2 Core Radius Error 


Fitted 


km 


km 

1-6 

0.0731 

3458 ±112 

- 22 

1-7 

0.0617 

3426 ± 

80 

- 54 

1-8 

0.0665 

3483 ± 

69 

+ 3 

1-9 

0.0765 

3538 ± 

63 

+ 57 

1-10 

0.0778 

3573 ± 

55 

+ 93 

1-11 

0.0763 

3548 ± 

47 

+ 68 

1-12 

0.2162 

3648 ± 

71 

+168 

1-13 

0.4410 

3766 ± 

93 

+286 

Table Fib: As Table Fla, but = 3. 


Degrees 

Fitted 

q 2 

Core Radius 
km 

Error 

km 

3-6 

0.0561 

3460 ±183 

- 20 

3-7 

0.0404 

3408 ± 108 

- 72 

3-8 

0.0530 

3516 ± 97 

+ 36 

3-9 

0.0620 

3598 ± 87 

+118 

3-10 

0.0586 

3638 ± 68 

+158 

3-11 

0.0664 

3587 ± 60 

+107 

3-12 

0.2116 

3723 ± 94 

+243 

3-13 

0.4230 

3875 ± 33 

+395 
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